Draft version February 2, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. ^jllj^A, 



in 
o 
o 

^3 



> 

m 
O 

\o 
o 
in 
o 

^: 
6 



X 



DISCOVERY OF X-RAY EMISSION FROM SUPERNOVA 1970G WITH Chandra: 
FILLING THE VOID BETWEEN SUPERNOVAE AND SUPERNOVA REMNANTS 

Stefan Immler^ and K. D. Kuntz^ 

Exploration of the Universe Division, X-Ray Astrophysics Laboratory, Code 662, 
NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
Draft version February 2, 2008 

ABSTRACT 

We report on the discovery of X-ray emission from SN 1970G in MlOl, 35 years after its outburst, 
using deep X-ray imaging with the Chandra X-ray observatory. The Chandra ACIS spectrum shows 
that the emission is soft 2 keV) and characteristic for the reverse shock region. The X-ray lu- 
minosity, £o.3-2 = (1-1 ± 0.2) X 10^^ ergs s""'^, is hkely caused by the interaction of the supernova 
(SN) shock with dense circumstellar matter. If the material was deposited by the stellar wind from 
the progenitor, a mass-loss rate of M = (2.6 ± 0.4) x 10~^ Mq yr~^ (uw/10 km s""'^) is inferred. 
Utilizing the high-resolution Chandra ACIS data of SN 1970G and its environment, we reconstruct 
the X-ray lightcurve from previous ROSAT HRI, PSPC, and XMM-Newton EPIC observations, and 
find a best-fit linear rate of decline of L oc with index s = 1.7 ± 0.6 over a period of 12-35 years 
after the outburst. As the oldest SN detected in X-rays, SN 1970G allows, for the first time, direct 
observation of the transition from a SN to its supernova remnant (SNR) phase. 

Subject headings: stars: supernovae: individual (SN 1970G) — stars: circumstellar matter — galaxies: 
individual (MlOl, NGC 5457) — ISM: supernova remnants — X-rays: general — 
X-rays: individual (SN 1970G, MlOl, NGC 4457) — X-rays: ISM 



1. INTRODUCTION 

Since the launch of the Chandra and XMM-Newton 
X-ray observatories, the number of supernovae (SNe) 
detected in their near aftermath has more than dou- 
bled (see Immler & Lewin 2003 for a review article)^. 
The high-quality X-ray spectra have confirmed the va- 
lidity of the circumstellar interaction models (see Frans- 
son, Lundqvist & Chevalier 1996 and references therein) 
which predict a hard spectral component (^ 10 keV) 
for the forward shock emission during the early epoch 
(^ 100 days) and a soft thermal component (^ 1 keV) 
for the reverse shock emission after the expanding shell 
has become optically thin. The soft emission component 
dominates the X-ray output of the interaction regions 
due to its higher emission measure and higher electron 
densities. This expected 'softening' of the X-ray spec- 
trum has been observationally confirmed for a number 
of young SNe, such as SNe 1978K (Schlegel et al. 2004), 
1979C (Immler et al. 2005), 1993 J (Zimmermann & As- 
chenbach 2003), 1999em and 1998S (Pooley et al. 2002). 

Where sufhcient data are available. X-ray lightcurves 
could be established over significant time scales (e.g., « 
25 years for SNe 1978K and 1979C; Schlegel et al. 2004, 
Immler et al. 2005; « 8 years for 1993J; Immler et al. 
2001, Zimmermann & Aschenbach 2003). However, little 
is known about the transition from a young SN into its 
supernova remnant (SNR) phase as there have been only 
few detections of intermediate-age (25-300 yrs) SNe in 
the radio and none in X-rays. While the X-ray emission 
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of SNe is dominated by the interaction of the shock with 
the ambient circumstellar matter (CSM), likely deposited 
by the progenitor's stellar wind, the emission from a SNR 
is thought to originate in the shocked interstellar medium 
(ISM). At an age of w 35 years, SN 1970G is the oldest 
SN detected in X-rays and closes this gap which allows 
one, for the first time, to witness the transition from a 
SN to its SNR phase. 

In this paper we report on the Chandra X-ray obser- 
vation of SN 1970G, as well as previous X-ray data from 
ROSAT and XMM-Newton. In S0 we briefly describe 
the data and analysis thereof and report on the X-ray 
spectrum and multi-mission X-ray lightcurve in ^ We 
discuss the results in the context of the CSM interaction 
model and put SN 1970G in the context of the evolution 
from a SN to a SNR in 21 followed by a summary in ^ 

2. DATA PROCESSING AND ANALYSIS 

SN 1970G has been observed with the Chandra ob- 
servatory as part of the 1 Ms observation of the host 
galaxy MlOl (NGC 5457; Kuntz 2005). Five individual 
Chandra ACIS observations (sequence numbers 600389, 
600389, 600389, 600389, 600389) from July 5-11, 2004, 
were merged into a single observation resulting in a 
cleaned and exposure-corrected on-source exposure time 
of 139.5 ks for the ACIS-S2 chip. Details of the data 
calibration and analysis are explained in Kuntz (2005). 
The high-resolution (w 0'.'5 FWHM) Chandra data were 
used to search for X-ray emission from the position of 
SN 1970G, and to examine the contamination by nearby 
X-ray sources. A Chandra ACIS-S2 image of SN 1970G 
and its environment is given in Fig. ^ 

We further constructed the long-term X-ray lightcurve 
of SN 1970G by extracting source counts from previ- 
ous ROSAT HRI (May 14 - Nov 23, 1996; sequence 
numbers RH600820N00, RH600820A01; merged expo- 
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Fig. 1.— Full (0.3-8 keV) band Chandra ACIS-S2 image of the 
region (size 1' X 1') around the position of SN 1970G. The image is 
smoothed with a Gaussian of 2'.' 5 (FWHM) and is given in logarith- 
mic gray-scale ranging from - 1.5 cts pix~^ (1 pix corresponds 
to 0'.'5). Contour Unes are 0.2, 0.4 (black), 0.6, 0.8, 1, 1.2, 1.4, and 
1.6 cts pix~^ (white). The radio position of SN 1970G (Cowan, 
Goss & Sramek 1991) north-east of the HII region NGC 5455 is 
marked by a white cross. 



sure time 176.1 ks) and ROSAT PSPC observations 
(Jun 8-9, 1991; sequence number RP600108N00; expo- 
sure 34.5 ks). Archival XMM-Newton EPIC observa- 
tions from June 4, 2002 (obs-id 0104260L01, 43.3 ks) 
July 23, 2004 (obs-id 0164560^01, 41.8 ks) January 8, 
2005 (obs-id 0212480201, 32.4 ks) were also analyzed ac- 
cording to standard analysis procedures described in the 
XMM-Newton ABC Guide version 2.01*. 

A short (2.8 ks) Einstein HRI observation 154 days af- 
ter the outburst of SN 1970G did not result in the detec- 
tion of the SN and is not included in this study. Archival 
ASCA observations were not used due to the large aper- 
ture (w 1' FWHM) and contamination by nearby X-ray 
sources. The ROSAT and XMM-Newton count rates 
were further compared to the published results by Wang, 
Immler & Pietsch (1999; ROSAT) and Jenkins et al. 
(2005; XMM-Newton). 

3. RESULTS 

A point-like X-ray source is detected at the radio 
position of SN 1970G (Cowan, Goss & Sramek 1991), 
with an offset (< 0'.'5) similar to the spatial resolu- 
tion and astrometry of the Chandra observation (see 
Fig. SN 1970G is spatially resolved and well sepa- 
rated from a nearby {~ 5") HII region (NGC 5455). Ex- 
posure, background, and aperture corrected counts ex- 
tracted from the position of SN 1970G gives a count rate 
of (4.3 ± 0.6) X 10"* cts s~^ While the hmited photon 
statistics does not allow a detailed characterization of 
the energy distribution of the recorded photons, spectral 
analysis shows that the emission is soft, with the bulk of 

* available from |http: / /heasarc.gsfc.nasa.gov/docs/xmm /abc/l 



the emission being confined to the ^ 2 keV band. 

Adopting a kT = 0.6 keV thermal plasma spectrum 
typical for the late-time (^ 100 days) emission of SNe 
and assuming a Galactic foreground column density with 
no intrinsic absorption (A^h — 1.16x10'^° cm^^; Dickey & 
Lockman 1990) gives a 0.3-2 keV flux and luminosity of 
/0.3-2 = (1.8 ±0.3) X 10-15 ergs cm-^ g-i and Lo.3-2 = 
(1.1 ±0.2) X 10'^^ ergs s^^, respectively, for a distance of 
7.2 Mpc (Stetson et al. 1998). 

Since no information is available about the spectral 
evolution of SN 1970G or other SNe at this age, we 
used the same spectral template to convert ROSAT and 
XMM-Newton count rates into fluxes and luminosities to 
construct the long-term (12-35 years) X-ray lightcurvc. 

An X-ray source is visible in the ROSAT HRI im- 
ages of MlOl at the position of SN 1970G (see Fig. 2 
in Wang, Immler & Pietsch 1999). The source, how- 
ever, is not included in the ROSAT HRI catalog since 
it is slightly below the employed detection threshold 
(S/N = 3.5; Wang, Immler & Pietsch 1999). Extrac- 
tion of source counts within the HRI aperture gives 
a 0.1-2.4 keV count rate of (2.9 ± 0.6) cts s'^, corre- 
sponding to a 0.3-2 keV band flux and luminosity of 
/0.3-2 = (6.7 ± 1.4) X 10-15 ergs cm'^ s^^ and Lo.3-2 = 
(4.2 ± 0.9) X 10^'' ergs s~i, respectively. 

An X-ray source is also listed in the ROSAT PSPC 
catalog at the position of SN 1970G (source P12; Wang, 
Immler & Pietsch 1999), with a signal-to-noise ratio of 
S/N = 4.5. Given that no other X-ray source is visible in 
the high-resolution Chandra image in the environment of 
SN 1970G within the PSPC aperture (« 25" FWHM), 
except for low-level emission of a neighboring HII region 
(NGC 5455, see Fig. P), it is safe to assume that the 
PSPC source is due to the SN itself. Assuming the emis- 
sion arises from SN 1970G, a (0.5-2 keV band) count 
rate of (1.3 ±0.1) x 10"*^ cts is inferred, correspond- 
ing to /o 3-2 = (1-3 ± 0.1) X 10"!'' ergs cm^^ s^^ and 
Lo.3-2 = (8.4 ± 0.9) X 1037 ergs s'^. 

A faint X-ray source is visible in each of the nine 
XMM-Newton EPIC images at low energies (< 2 keV). 
No photons were recorded in the hard (2-6 keV) band. 
The source, however, is extended (« 20") and shows no 
concentration of the recorded photons at the radio po- 
sition of SN 1970G. Since it is likely that the emission 
is contaminated by the nearby HII region and other soft 
emission components inside the galactic disk, we only 
use the count rates, fluxes, and luminosities as upper 
limits. Discrepancies between the published PN and 
MOS count rates for source #44 (factor xl6.5; Jenkins 
et al. 2005) for one of the pubhshed XMM-Newton data 
(obs-id 0104260101) makes a comparison problematic. 

We calculated the mass-loss rate of the progenitor as 
a function of the stellar wind age using the relationship 
= 4/(7rr7i2)A(T) x {M/v^)^ x (wst)~\ where m is 
the mean mass per particle (2.1 x 10^^ kg for a H+He 
plasma), A(T) the cooling function of the heated plasma 
at temperature T, M the mass-loss rate of the progenitor, 
Vyj the speed of the stellar wind blown off by the progen- 
itor, and Vs the speed of the outgoing shock (see Imm- 
ler, Wilson & Terashima 2002). A constant shock veloc- 
ity of Vs = 9, 000 km s^^ similar to other core-collapse 
SNe, such as SN 1979C (Marcaide et al. 2002) was as- 
sumed. An effective (0.3-2 keV band) cooling function 
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TABLE 1 
X-Ray Properties of SN 1970G 



Epoch 


Instrument 


Count Rate 


/0.3-2 


Lo.3-2 


M 


[years] 




[10-4 cts s-i] 


[10—-*-^ ergs cm"^ -"-] 


[10^'' ergs s-i] 


[10-4 yr- 


11.9 


ROSAT PSPC 


13.0 ±1.3 


13.4 ± 1.4 


8.4 ±0.9 


4.2 ±0.4 


16.8 


ROSAT HRI 


2.9 ±0.6 


6.7± 1.4 


4.2 ±0.9 


3.4 ±0.7 


31.9 


XMM-Newton EPIC 


< 48.3 


< 10.2 


< 6.3 


< 5.8 


34.0 


XMM-Newton EPIC 


< 80.6 


< 17.0 


< 10.6 


< 7.8 


34.5 


XMM-Newton EPIC 


< 119.3 


< 25.2 


< 15.6 


< 9.4 


34.7 


Chandra ACIS 


4.3 ±0.6 


1.8 ±0.3 


1.1 ± 0.2 


2.6 ±0.4 



of A = 3 X 10^^^ ergs cm'^ s""'^ for an optically thin ther- 
mal plasma with a temperature of 10^ K was adopted, 
which corresponds to the assumed thermal plasma tem- 
perature. Assuming different plasma temperatures in 
the range 0.5-1 keV would lead to changes in the emis- 
sion measure of ^ 10%. Key observational properties of 
SN 1970G are listed in Tabled 

4. DISCUSSION 

Previous X-ray observations of SN 1970G with 
ROSAT, ASCA, and XMM-Newton lacked the spatial 
resolution needed to separate the SN from the nearby 
(« 5") HII region NGC 5455. Using the sub-arcsec 
Chandra imaging capability, in combination with the 
deep exposure time (140 ks), we detect soft 2 keV) 
X-ray emission from SN 1970G with a luminosity of 
Lq.3-2 ~ 1 X 10"^^ ergs s"^, 35 years after its outburst. 
The source is spatially resolved from the nearby HII re- 
gion NGC 5455 (« 5"), which is observed at a <'20% flux 
level compared to SN 1970G. Contamination of the previ- 
ous ROSAT data of SN 1970G with the HII region should 
therefore be minimal, especially when kept in mind that 
at earlier epochs SN 1970G was significantly brighter, 
while the HII region presumably has a constant flux. 

Given the information from the high-resolution spatial 
imaging with Chandra, we recover X-rays from SN 1970G 
in the ROSAT HRI and PSPC observations. SN 1970G 
is not conclusively detected in three recent (2002-2005) 
XMM-Newton observations due to the extent of the 
source (w 20") as imaged by the PN and MOS instru- 
ments, and positional offsets larger than the instrumental 
uncertainties (« 5"). 

In combination, an X-ray rate of decline of L oc 
with index s ~ 1.7 ± 0.6 is similar to other X-ray SNe 
(e.g., SN 1999em, Pooley et al. 2002; SN 1986J, Houck 
et al. 1998). This strengthens the case that the emission 
is due to the SN itself, even though a variable source in 
the environment of NGC 5455 cannot be entirely ruled 
out. The lack of X-rays above 2 keV confirms that the 
emission is soft and of possible thermal original. Such a 
soft emission is expected for the late emission of a SN, 
where the shocked emission originates from the reverse 
shock region. 

Assuming that the emission arises from shock 
heated plasma deposited by the progenitor's stel- 
lar wind, a mass-loss rate of M — (2.6 ±0.4) x 
10^^ Mq yr^-^ (ww/10 km s^^) is inferred for the 
late epoch observed with Chandra. A similar mass- 
loss rate was inferred from radio observations [2 x 
10-5 Mo yr-i (ww/10 km s^^); Weiler 1993]. The 



mass-loss rates for the earlier epoch are slightly higher 
and might indicate an evolution in the mass-loss rate 
history of the progenitor (see Tab. It should be 
noted, however, that the lack of knowledge about the 
shock front velocity precludes a more detailed study. If 
the shock front experienced a significant deceleration, 
as expected when the shock has swept up a significant 
amount of GSM, the lower inferred mass-loss rates for 
the later epochs would be even lower (see the discus- 
sion about the mass-loss rate history of SN 1979C in 
Immler et al. 2005). However, given the uncertainties in 
the derived X-ray luminosities due to a possible resid- 
ual contamination of the ROSAT and XMM-Newton 
data caused by their larger point-spread functions (PSF; 
« 5" FWHM for the ROSAT HRI; 5"~6" FWHM for 
the XMM-Newton EPIC instrument; « 25" FWHM for 
the ROSAT PSPC), the apparent increase of the mass- 
loss rate at earlier epochs is not regarded as significant. 
Interestingly, the mass-loss rates are higher for the in- 
struments with the largest PSF (EPIC and PSPC). If 
the mass-loss rates inferred from the instruments with 
the smallest PSF are compared [M — (3.4 ± 0.7) and 
(2.6 ± 0.4) X 10-5 Mq yr-i (?;w/10 km 3-^) for the 
HRI and ACIS, respectively], no statistically significant 
differences are observed. We therefore conclude that 
SN 1970G shows a mass-loss rate which is consistent with 
being constant over a period of 12-35 years after the 
peak optical brightness (July 30, 1970; Detre 1970), cor- 
responding to « 11,000-31,000 years in the stellar wind 
age using the conversion t„ ~ ivs/v^ and assuming a 
stellar wind speed of 10 km s- and a shock front veloc- 
ity of 9,000 km s-\ 

The measured mass-loss rate for SN 1970G is simi- 
lar to those inferred for other Type II SNe, which typi- 
cally range from 10-^ to 10-** Mq yr~^ (see Immlcr & 
Lewin 2003). This is indicative that the X-ray emission 
arises from shock-heated CSM deposited by the progeni- 
tor rather than shock-heated ISM, even at this late epoch 
after the outburst. 

In order to explain the X-ray luminosities of SNRs such 
as Gas A (Lo.3-2 = 2.8 x 10'^'' ergs s-^ at an age of 
« 320 years), mass-loss rates in excess of 10-** M© yr-^ 
would be required in the shocked stellar wind scenario. 
The high required amounts of deposited material could 
not be provided by a sustained stellar wind since the 
wind would interact with the ambient ISM at these large 
radii (^ 10^^ cm) from the site of the explosion. Instead, 
the X-ray emission from SNRs is thought to originate 
in radiative cooling of the shock-heated ISM which has 
substantially larger densities. The relatively low inferred 
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Fig. 2.— Soft (0.3-2 keV) band X-ray luminosities of all SNe 
detected to date (filled squares) and historical SNRs (filled trian- 
gles) as a function of age (in units of years). The X-ray lightcurve 
of SN 1970G is marked by filled diamonds with error bars (from 
left to right: ROSAT PSPC, ROSAT HRI and Chandra ACIS). 
XMM-JVewton EPIC upper limits are indicated by arrows. 



mass-loss rate required to produce the observed X-ray 
luminosity of SN 1970G indicates that the SN shock has 
not yet reached the ISM. Consequently, the SN 1970G 
has not yet reached its SNR phase. 

We therefore conclude that at an age of 35 years, 
SN 1970G is the oldest SN observed in X-rays, and pro- 
vides important constraint as to when the transition from 
a SN to its SNR phase occurs, which is likely at a later 
stage (« 50-100 years) after its outburst. 

In Fig. 13 we plot the X-ray evolution of all SNe de- 
tected to date^ as well as the X-ray luminosities of histor- 
ical Galactic, LMC, SMC (as observed with Chandra^), 
M31 (Supper et al. 2001; Kong et al. 2003), and M33 
(Ghavamian et al. 2005) SNRs. A comparison of the X- 
ray luminosities shows that the transition from a SN to 
a SNR appears to be rather smooth. If the evolution of 
SN 1970G continues as observed to date, it will be simi- 
lar to X-ray faint Galactic SNRs such as SNRs 1181 and 
1006. X-ray bright SNe such as SNe 1978K and 1979C, 
which show no sign for an X-ray evolution out to large 
distances (« 10^^ cm) at an age of « 25 years after their 
outburst, can evolve into luminous SNRs such as Cas A 

^ http: / /Iheawww. gsfc.nasa.gov/users/immler/supernovae_list. html 
® http:/ /www. astro. psu.edu/users/green/Main/main5. html 



due to their strong interaction with their dense ambient 
GSM, or possibly attain the X-ray luminosities of the 
brightest SNRs observed (such as W49B or N132D) in 
case of a dense ISM interaction. 

Clearly, X-ray detections and monitoring of 
intermediate-age (around 100 years) SNe are needed to 
address these questions and to study the transition from 
a SN into a SNR in detail. Deep galaxy surveys with 
Chandra and XMM-Newton and future X-ray missions 
such as Const-X and XEUS provide the capabilities to 
detect these intermediate- age SNe in the near future 
and to fill the void in the evolution of a SN to a SNR. 

5. SUMMARY 

Based on our deep (140 ks) X-ray imaging of MlOl 
with the Chandra observatory we present the detection 
of SN 1970G in X-rays. Key results are: 

• SN 1970G is detected with a luminosity of 
io.3-2 = (1-1 ± 0.2) X 10^^ ergs s-\ 35 years after 
its outburst in our deep Chandra observation of MlOl. 

• SN 1970G is recovered in previous ROSAT HRI and 
PSPG observations and shows a best-fit linear rate of 
decline of L (x with index s = 1.7 ± 0.6 over the 
observed period of 12-35 yrs after its outburst, consistent 
with recent XMM-Newton observations. 

• The emission of SN 1970G is soft {<, 2 keV) and indi- 
cates that it originates in the reverse shocked region, as 
expected for the late emission in the evolution of a SN. 

• A mass-loss rate of M = (2.6 ± 0.4) x lO"'^ M© yr^^ is 
inferred consistent with being constant over a period of 
11,000-31,000 years in the evolution of the progenitor. 

• At an age of 35 years, the SN has not yet reached its 
SNR phase, which is likely to occur at a later stage in 
the evolution (w 50-100 years). 

• Comparison of the X-ray lightcurves of all SNe detected 
to date with the X-ray luminosities of SNRs indicates a 
smooth transition of SNe into their SNR phases. 

This paper is based on observations obtained with 
XMM-Newton, an ESA science mission with instruments 
and contributions directly funded by ESA Member States 
and NASA. KDK was supported by the Chandra grant 
SAO GO-5600587. The authors wish to thank R. Petre 
for helpful suggestions and comments. 
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